Pressure spectrum and structure function in homogeneous turbulence 
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The pressure spectrum and structure function in homogeneous steady turbulence of an incom- 
pressible fluid is studied using direct numerical simulation. The resolution of the simulation is up to 
1024 3 and the Taylor microscale Reynolds number R\ is between 38 and 478. The energy spectrum 
is found to have a small but finite inertial range followed by a bump at large wavenumbers. The 
Kolmogorov constant K is found to be 1.66 ± 0.08. The pressure spectrum also has a small but 
finite inertial range of P(k) = _B p e 4,/3 fc _7/ ' 3 followed by a bump of nearly k~ 5 ^ 3 range at higher 
wavenumbers. Both scaling ranges match at a crossover wavenumber, k p , which is a characteris- 
tic wavenumber for the pressure gradient. The constant B p is found to be about 8.34 ± 0.15 for 
R\ — 460. Its non-universality is discussed. The second order pressure structure function, computed 
in terms of the fourth order velocity structure functions, agrees well with that obtained by the direct 
measurement over separations ranging between the inertial and dissipation scales. 
PACS numbers: 47.27.Ak, 47.27.Jv, 47.27.Gs, 05.20.Jj 



The pressure spectrum in a turbulent incompressible 
flow is defined as (p 2 ) = J °° P(k)dk. Kolmogorov's the- 
ory predicts that 



P{k) 



L 1 < k < rf 



(1) 



where <fr(x) is a non-dimensional function, v is the kine- 
matic viscosity, e is the average rate of energy dissipation 
per unit mass, L is the integral scale of turbulence, rj is 
the Kolmogorov scale, and B p is a non-dimensional con- 
stant of order one. The fluid density is assumed to be 
unity throughout this paper. 

There have been many studies of the pressure spec- 
trum. [i]-[2C|l Some of the experiments have shown 
that P(k) oc fc~ 7 / 3 , H or equivalently D p — 
{{p(x + r) — p(a?)) 2 ) oc r^. ]12|,f20| Others have re- 
ported that r 4 / 3 is not observed. jf3| Recent DNS's with 
large scale forcing have found that the pressure spectrum 
is approximately proportional to fc -5 / 3 , unlike k~ 7 ^ 3 , in 
the wavenumber range where the energy spectrum scales 
close to fc~ 5 / 3 . |^fj|-|l8[| Gotoh and Rogallo conjectured 
that the observed k~ 5 ^ 3 scaling for P(k) is a bump, as 
observed for the energy spectrum, and that P(k) scales 
as fc~ 7 / 3 in the lower wavenumber range. This implies 
that a wider inertial range, compared to the energy spec- 
trum, is necessary for K41 scaling of P(k). There 
seems to be no agreement about the scaling of the pres- 
sure spectrum when compared to the case of the energy 
spectrum. 

We have performed a series of DNS's of incompress- 
ible homogeneous isotropic turbulence using a resolution 
of up to N = 1024 3 . The DNS was designed to gen- 
erate a wider inertial range and higher Reynolds num- 
bers. The range of the Taylor microscale Reynolds num- 
ber R\ = u\jv is between 38 and 478, where u is the 



root mean square of turbulent velocity and A is the Tay- 
lor microscale. The characteristic parameters of the DNS 
arc listed in Table I. The code uses the pseudo Fourier 
spectrum and 4th order Runge Kutta Gill methods. Ran- 
dom forcing, Gaussian and white in time, is applied to 
the lower wavenumbers. A statistically steady state was 
confirmed by observing the time evolution of the total 
energy, the total enstrophy and the skewness of the longi- 
tudinal velocity derivative. The statistical averages were 
taken as the time average over tens of turnover times for 
lower Reynolds numbers and over a few turnover times 
for the higher Reynolds numbers. The data of the highest 
Reynolds number, R\ = 478, were obtained as short time 
average (about 0.34 eddy turnover times) during the pas- 
sage toward steady state rather than over a statistically 
steady state. The resolution condition k max rj > 1 is sat- 
isfied for most runs, but that of the case when R\ = 460 
is slightly less than unity. We believe that this does 
not adversely affect the energy and pressure spectra re- 
sults in the inertial range. Computations with R\ < 284 
were performed using a Fujitsu VPP700E vector parallel 
machine with 16 processors at RIKEN. Simulations us- 
ing higher R\ were performed on a Fujitsu VPP5000/56 
with 32 processors at the Nagoya University Computa- 
tion Center. 

The energy spectra in Kolmogorov units (multiplied by 
(/cr/) 5 / 3 ) are shown in Fig.|l|. Collapse of curves at vari- 
ous Reynolds numbers is very satisfactory, although the 
curves with R\ > 284 have appreciable rise of E(k) near 
the high wavenumber boundary. The Kolmogorov con- 
stant K was measured in the range of 0.008 < kr\ < 0.04 
in which the average energy transfer flux function -ZT(fc) /e 
is nearly flat and close to unity (figure not shown). The 
value of K 



K = 1.66 ±0.08, 



(2) 
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is very close to the value of 1.62 obtained in previous 
experiments and DNS's, [^l],|2| and to the value of 1.72 
obtained using the Lagrangian spectral theory (LRA). 
p3| , p4| There is a small spectral bump at wavenumbers 
near krj w 0.2, as observed in other DNS's. |22| 

Figure 2 shows the pressure spectra in terms of the K41 
scaling, Eq.JIJ], multiplied by (/cry) 7 / 3 for various Reynolds 
numbers. Fig. ^ is a close-up of the curves for higher 
Reynolds numbers. For R x < 300 there is no plateau 
in the curves. However, for R x larger than 400, there 
appears to be a small plateau of finite length. There is 
a bump, with a peak value of about 17, in P(k) near 
krj = 0.2 which is more appreciable than that in the en- 
ergy spectrum. The left part of the bump consists of 
a finite ramp. For R\ — 284, the slope of the ramp is 
close to 2/3, indicating that P(k) oc fc~ 5 / 3 . The slope 
gradually decreases as the Reynolds number increases. 
It is this part which the previous DNS's have observed 
as P{k) oc fc~ 5 / 3 . 

The curves obtained for R x = 387, 460 and 478 indi- 
cate that P(k) approaches the fc _7//3 spectrum over the 
range of 0.008 < k-q < 0.04. The value B p is 

Bp = 8.34 ±0.15, (3) 

shown in Fig. 3 as a horizontal line. Taking into account 
the relatively short length of the averaging time, the er- 
ror bar for B p is a few times larger than 0.15. Pullin 
obtained B p = 1.325K 2 using the joint Gaussian hypoth- 
esis for the 4th order velocity structure functions. With 
K = 1.66, Bp = 3.65, which is smaller than the present 
DNS value. This is consistent with the fact that P(k) 
is larger than Pa(k), computed from the Gaussian ran- 
dom velocity field with the same energy spectrum as that 
of the actual turbulence field. |]]|,[l8|] P urmr suggested 
B p ks 7 using DNS data with N = 128 3 at R x = 77.5. 
Pullin estimated that B p »s 2.14 — 3.65 using Lundgren's 
stretched spiral vortex model. 

The collapse of the P(k) curves for all R x s is not 
as good as the energy spectrum, even in the dissipa- 
tion range. The collapse of the pressure spectra is im- 
proved when the normalized pressure gradient variance, 
Fv P = ((Vp) 2 )^ 3 / 2 ;/ 1 / 2 , is included in the scaling for 
P{k): 

P(fc) = F Vp e 3 /V/V l (M» ( 4 ) 

where <fti(x) is a non-dimensional function. |]l9| , p5| Fig- 
ure ^ presents P(k) using Eq.[Q, and clearly shows that 
the scaling of P(k) in the high wavenumber range is bet- 
ter than the scaling using Eq.jl]]. The inset shows the 
variation of F^ p against the Reynolds number. F^ p is 
a monotonically increasing function of R x that becomes 
very weakly dependent on R x as R x becomes large. It 
should be noted that although the insensitivity of F-y P to 
R x is consistent with Batchelor's Gaussian theory for the 
pressure, its value is considerably larger than the value 



corresponding to the Gaussian theory, F^ p . Jl],|l8| This 
insensitivity of Fy p at large R x implies that the collapse 
of P(k) with R x > 284 is little affected by Fy p . However, 
there still remains a weak Reynolds number dependence 
of P(k) in the inertial range, causing the pressure spec- 
trum to be non-universal. Close inspection of P(k) in 
the fc -7 / 3 range shows that the factor F\j p of Eq.Jpl im- 
proves the collapse of the curves (Figs. || and In 
this range, tj>i{krj) — C p (kri)~ 7 / 3 , where C p is a non- 
dimensional constant of the order one. The constant B p 
is related as B p = Fy p C p , so that B p becomes weakly 
dependent on Reynolds number, while C p is not. It is 
reasonable, in this sense, to regard C p as a more univer- 
sal constant than B p . Using the values of B p and F\/ p 
we obtain 

C p = 0.707 ±0.1. (5) 

The non-universality enters the pressure spectrum 
through F\j p {R x ) as a function of R x . Therefore, Fy p is a 
key parameter for the second order statistics of pressure. 
The Reynolds number dependence of i*vp is attributed 
to the coherent structure of the source term field in the 
Poisson equation for the pressure. [|l8|]25| l 

Transition to the nearly k~ 5 ^ 3 range occurs at kr\ « 
0.03 for R x > 387, which corresponds to k p X p « 1.5, 1.8 
and 1.6 for R x — 387, 460 and 478, respectively. Here, A p 
is a characteristic length scale for the pressure gradient 
defined by ((dp/dx) 2 ) = u 4 /Xp. This is analogous to the 
Taylor microscale. Thus the crossover scale between the 
/c~ 7 / 3 and the nearly A: -5 / 3 ranges is about k p — A" 1 . 

The pressure spectrum at low wavenumbers scales as 

P(k) = ^ 3 L 7 / 3 (j> 2 {kL), (6) 

where L is the integral scale. The curves of P(k) at this 
range collapse reasonably well into one curve (figure not 
shown). The scaling of Eq.Q matches Eq.j(| in the k~ 7 1 3 
range. 

Hill and Wilczak have derived an expression for the 
pressure structure function D p (r) in terms of the fourth 
order velocity structure functions, assuming a homoge- 
nous and isotropic velocity field: 

1 4 f°° 

D p (r) = --L(r) + -r 2 J y- 3 [L(y) + T(y) - &M(y)]dy 

where L(r) = (5u(r) 4 ),T(r) = (Sv(r) 4 ), M(r) = 
(8u{r) 2 8v{r) 2 ^ , and 5u(r) and Sv(r) are the longitudi- 
nal and transversal velocity differences, respectively. |Tq1 
Isotropy has been examined in terms of the kinematic 
conditions for the second and third order moments of the 
longitudinal and transversal velocity increments. The re- 
lations are well satisfied for both R x = 387 and R x = 460. 
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(The relative error in the equation for the third order mo- 
ments is less than about 10% at r/rj = 200.) 

Figure || shows plots of L(r),T(r) and M(r) at R\ = 
387 and 460. There is a straight line contained along 
each curve between 50 < r/rj < 400. The straight line 
showes the slope of 1.28, the value predicted by She and 



Leveque. 27 An examination of the local exponent of 
the curves showed that the inertial range exponents de- 
fined by L(r) ex r^ 4 , T(r) oc r^ 4 and M(r) oc r^ 4 have 
plateau between 90 < r/77 < 200 and 1.24 < cj < < 
Cf < 1.32 for R x = 460. §(| The detailed analysis will 
be reported elsewhere. J2q 



D p (r) computed by Eq.jTj is compared with values ob- 
tained using direct measurement in Fig. ^. The curves 
for Rx = 460 are shifted upward by one unit for clarity. 
Agreement of the curves for r/rj < 150 is satisfactory. It 
is reasonable to expect that the isotropy of the fourth 
order velocity structure functions in the limit of r — > 0, 
the fundamental assumption for the derivation of 0, is 
well satisfied. The pressure gradient variance and the 
pressure structure function at small separations can be 
examined in terms of Eq. j?J . |h],^(| In the inset of Fig^, 
the values of FVp computed using D p (r)r~ 2 in the limit 
of r — > are found to be very close to the value of Fy p 
obtained by direct measurement. 

The straight line in FigJ| indicates the r 1 slope. The 
slope of D p (r) is very close to unity and between 2/3 and 
4/3. The scaling of D p (r) requires a much longer scale 
separation than the case using wavenumber space. 
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is the length of time average. 
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Fig.l Scaled energy spectra, e 
K = 1.66 ±0.08. 



-1/4 



v- 5 ' 4 {kr,f' 3 E{k). 



Fig. 4 Scaling of the pressure spectra with the fac- 
tor F Vp , F^e- 3/4 u- 7/4 (fcr;) 7/3 P(fc). The lines are 
the same as those in Fig. 2. The inset is the vari- 
ation of -FVp against R\. Squares are F\? p computed 
by lim, — >o D p (r)/r 2 , and other symbols are experimental 
data by Voth et al. pi 



Fig. 2 K41 scaling for the pressure spectra, e 
(/b?) 7/3 P(fc). 



--3/4-7/4 



Fig.3 Close up of e" 3/4 u~ 7/4 {krj) 7/3 P(k) for higher 
Reynolds numbers. A short straight line shows the slope 
of k 2 / 3 and a horizontal line indicates B„ = 8.34. 



Fig. 5 Fourth order structure functions of velocity in- 
crements. R\ = 387 and R\ — 460 (thick lines). A 
straight line shows the slope of 1.28. 



Fig.6 Comparison of D p {r) with UkW.^R x = 387 
and R\ — 460 (thick lines). A straight line shows the 
slope of r 1 . The curves for R\ — 460 are shifted upward 
by one unit for clarity. 
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